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We study the low-temperature isothermal magnetic hysteresis of cubical and spherical nanoparti¬ 
cles with ferromagnetic (FM) core - antiferromagnetic (AF) shell morphology, in order to elucidate 
the sensitivity of the exchange bias effect to the shape of the particles and the structural imper¬ 
fections at the core-shell interface. We model the magnetic structure using a classical Heisenberg 
Hamiltonian with uniaxial anisotropy and simulate the hysteresis loop using the Metropolis Monte 
Carlo algorithm. For nanoparticles with geometrically sharp interfaces, we find that cubes exhibit 
higher coercivity and lower exchange bias field than spheres of the same size. With increasing inter¬ 
face roughness, the shape-dependence of the characteristic fields gradually decays and eventually, 
the distinction between cubical and spherical particles is lost for moderately rough interfaces. The 
sensitivity of the exchange bias field to the microstructural details of the interface is quantified by 
a scaling factor (6) relating the bias field to the net moment of the AF shell ( H e b = bMAF + H 0 ). 

Cubical particles exhibit lower sensitivity to the dispersed values of the net interfacial moment. 

PACS numbers: 75.50.Tt, 75.75.Fk, 75.75.Jn 
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I. INTRODUCTION 

Magnetic nanoparticles composed of two coupled mag¬ 
netic phases with distinct magnetic properties, such as 
magnetic anisotropy, are presently extensively studied, 
owing to their potentials for a wide variety of technolog¬ 
ical applications ranging from magnetic recordings to 
biomedicine £ In particular, nanoparticles with the most 
commonly achieved core-shell morphology are grown ei¬ 
ther by vacuum chamber condensation methods^ or more 
often by colloidal chemistry methods^ The latter meth¬ 
ods are low-cost and very versatile, producing nanoparti¬ 
cles with a very narrow size distribution and a variety of 
combined materials for the core and the shell. Owing to 
the flexibility offered by synthetic methods, nanoparticle 
shape has emerged as a promising tool for tailoring the 
magnetic properties at the nanoscale^ Prominent exam¬ 
ples include engineering of nanoparticle shape to obtain 
a better response for the biomedical imaging^ or mag¬ 
netic hyperthermia applications £ Nanoparticles with a 
ferromagnetic (FM)-core / antiferromagnetic (AF)-shell 
morphology are well known to exhibit the exchange bias 
effect after been field-cooled in a field^^. Two charac¬ 
teristic features of the effect being the horizontal shift of 
the hysteresis loop in the direction opposite to the cool¬ 
ing field and the widening of the loop, expressed by an 
increase of the coercivity relative to the bare FM parti¬ 
cle. In most cases, the exchange bias phenomenon arises 
from the pinning of magnetic moments at the interface 
between the two materials : 10 ’ 11 Thus the material proper¬ 


ties, the size and the atomic structure of the interface are 
crucial factors determining the existence and the strength 
of this phenomenon. 

Despite the abundance of experimental and theoreti¬ 
cal works focusing on the role of the core size and shell 
size in controlling the exchange bias phenomenon in bi- 
phased nanoparticlesiS, less attention has been payed 
so far to the interplay between exchange bias and par¬ 
ticle shape. Experimental works have reported on the 
exchange bias effect in core-shell nanoparticles with cu¬ 
bical shape^i^F— and more recently nanoparticles with 
more complicated geometrical shape have been synthe¬ 
sized, including octahedra, cuboctahedra and octopods^ 
trianglesdimers and flowers^, and their magnetic 
properties have been reported. At the same time, the¬ 
oretical works on the magnetic properties of core-shell 
nanoparticles with non-spherical shape are a few, 17,18 de¬ 
spite the necessity to unravel the interplay between ex¬ 
change bias and particle shape. In addition to the role 
of particle shape, the importance of the FM-AF inter¬ 
face quality in controlling the magnetic properties of the 
exchange-coupled bi-phased nanoparticle has been recog¬ 
nized. In a recent work, Mao et al^ demonstrated that 
the presence of non-magnetic atoms (vacancies) at the 
core-shell interface degrades the exchange bias effect in 
FM/AF nanoparticles. Contrary to this trend, Evans et 
al2& demonstrated that strong exchange bias can develop 
in weakly-coupled core-shell Co/CoO nanoparticles due 
to the presence of interface roughness. Finally, in recent 
experimental work by Juhin et i direct evidence on 


2 


TABLE I: Structural parameters of modeled nanoparti- 
cles with ideal interfaces. R = radius (sphere) or half¬ 
edge (cube) R c = core radius or core half-edge, A=total 
number of spins (atoms), N c = core spins, N c if=core inter¬ 
face spins, iV s / 1 i/=shell interface spins, lV lm c=uncompensated 
spins, S=sphere, and C=cube. 


Particle 

R/a 

R c /a 

N 

N c 

Ncif 

Nshif 

Nunc 

S710 

10.5 

7.5 

4945 

1189 

602 

762 

42 

S58 

8.0 

5.0 

2109 

257 

258 

410 

30 

C58 

8.0 

5.0 

4913 

729 

602 

866 

2 

C58r* 

8.0 

5.0 

4505 

549 

694 

642 

38 


(*) Cube faces along the (110), (110) and (001) planes. 


the existence of an interdiffused layer at the interface of 
core-shell nanoparticles was provided, thus showing the 
way to obtaining atomic scale information on the quality 
of the buried core-shell interface. 

Provided that the overall shape of the nanoparticle 
and the atomic scale quality of the interface are two 
crucial factors determining the strength of the exchange 
bias effect, in the present work, we implement numeri¬ 
cal simulations to study the interplay between these two 
factors in the case of cubical and spherical FM(core)- 
AFM(shcll) nanoparticles. Our main conclusions are that 
cubes show higher coercivity and lower exchange bias 
field than spheres, while even moderate interface rough¬ 
ness washes out the shape-dependence of the exchange 
bias effect. Furthermore, in assemblies of cubical par¬ 
ticles a wider dispersion of exchange bias field values is 
expected, due to the lower compensation of flat core-shell 
interfaces. 


II. ATOMISTIC MODEL 

The atomic sites of a nanoparticle are generated by 
cutting a sphere of radius R or a cube of edge L from an 
infinite simple cubic (SC) lattice with lattice constant a. 
The center of the cube or the sphere coincides with the 
central atomic site of the nanoparticle. For both particle 
shapes, the core has the same geometrical shape as the 
whole particle and the shell thickness is taken t s h = 3 a. 
The interface region is considered two atomic layer thick, 
one atomic layer on either side. The structural parame¬ 
ters of the nanoparticles studied in the present work are 
summarized in Table Hj The particular choice of sizes is 
made in order to emphasize the importance of shape ef¬ 
fects through comparison of cubes and spheres with the 
same size (558, (758, C58r ) or approximately the same 
volume, i.e. number of atoms (5710, (758, (758r). The 
atomic scale structure of the core-shell interface is not 
easily controlled experimentally. In the most common 
case of surface oxidized metallic nanoparticlesii, the in¬ 
terface profile is formed by a diffusion-like process and 
is expected to have a highly complex structure contain¬ 


ing various types of structural and substitutional defects. 
Here, we approximate the complex interface profile by a 
model of atomic scale roughness at the core-shell inter¬ 
face, which forms by the interpenetration of the magnetic 
phases. For simplicity we assume that roughness is con¬ 
fined within the two atomic layers of the interface region. 
Numerically, the rough interface profile of a spherical par¬ 
ticle is realized by sampling the local radius at each site 
of the core-interface region from a Gaussian distribution 
with mean value equal to the core radius ( R c ) and vari¬ 
ance a. Correspondingly, for cubical particles, the local 
edge size is Gaussian distributed. Typical particles with 
rough interfaces are drawn in Fig [T] This atomistic model 
of interface disorder minimizes the formation of isolated 
substitutional defect sites^H and preserves to a reasonable 
degree the spatial separation of the two magnetic phases. 



FIG. 1: (Color on line). Mid-plane of a cubical (C58) and 
a spherical (5710) particle with rough interfaces. Roughness 
dispersion is a — 0.5a in both cases. 

The magnetic structure is described by a classical 
Heisenberg Hamiltonian with first nearest neighbor ex¬ 
change interactions and uniaxial anisotropy along the 
^-axis, namely H = J ij s i ' s j - J2i K i S iz ~ 

where S, is a classical spin vector on site i. 
The exchange energy constant Rj takes the values Jfm, 
Jaf and J lnt depending on whether sites i, j belong both 
to the FM region, both to the AF region, or one to the 
FM and the other to the AF region, respectively. Simi¬ 
larly, the anisotropy energy constant Ki takes the values 
Kfm and I\af , depending on the type of material on site 
i. Magnetostatic (dipolar) energy terms are neglected in 
the previous Hamiltonian, because the size of nanoparti¬ 
cles studied here, is well below the single-domain limit.- 22 
In our simulations we use dimensionless energy parame¬ 
ters that have been previously introduced-^ to describe 
Co(core)/CoO(shell) nanoparticles. In particular, we de¬ 
fine Jfm as the unit of energy ( Jfm = 1) and then 
Jaf = —0.5 JFM,Jint = —0.5 Jfm, Kfm = 0.1 Jfm, 
and Kaf — 0.5 Jfm- In this set of parameters the AF 
anisotropy value seems to be too large ( Kaf ~ Jaf)^- 
However, we should note that since the AF shell is very 
thin (t s h = 3a), one expects the existence of strong sur¬ 
face anisotropy, which justifies the above choice of param¬ 
eters. In our dimensionless units, temperature T is mea¬ 
sured in units of Jfm/^-b and the magnetic field strength 
H in units of Jfm/ g^B, where gs is the Bohr magneton 
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and g the Lande factor. 

It is well known^iii that in order to observe shifted 
hysteresis loops, the composite nanoparticles have to be 
field cooled from a high temperature (T Tc) to a 
low temperature (T -C Tjv) in a relatively weak field 
(. H coo i < H sa t) prior to sweeping the applied magnetic 
field (—H coo i < H < +H coo i). At the end of the FC pro¬ 
cess the FM moments are aligned along the field direction 
while the net AF moment is frozen parallel ( J int > 0) 
or antiparallel (,/ mt < 0) to the cooling field produc¬ 
ing an exchange-bias field on the FM. Due to the uniax¬ 
ial symmetry of our model Hamiltonian and the strong 
anisotropy of the AF material, the exact spin configura¬ 
tion of the AF shell at the end of the FC process can 
be approximately described by an Ising-type AF con¬ 
figuration, which is either type-I (-1-1—) or type-II 

(—|-h). Adopting this approximation for the FC state 

of the AF, we avoid the time-consuming simulation of the 
actual FC process. Then the hysteresis loop simulation 
starts using as the initial spin configuration of the AF, 
the one (type-I or II) that minimizes the total energy of 
the particle under the applied H coo i fields. Tests runs 
have shown that this approximation to the FC state, pro¬ 
vides configurationally average loops almost identical to 
those obtained by the exact calculation of the FC pro¬ 
cess. Hysteresis loops are simulated using the standard 
Metropolis Monte Carlo algorithm with single spin up¬ 
dates and 10 4 initial Monte Carlo steps per spin (MCSS) 
for thermalization followed by 10 4 MCSS for thermal av¬ 
eraging. Sampling over thermal disorder is performed 
every r = 10 MCSS to minimize correlations between 
sampling points. The percentage of accepted spin moves 
is kept close to 50% by adjusting the width of the spin 
moves. A field sweep rate A H/MCSS = 10 -5 is used, 
and is kept constant throughout the present work, in or¬ 
der to point out changes in the loop characteristics aris¬ 
ing solely by the structural and morphological properties 
of the particles and not the sampling time. For a given 
particle shape, size and dispersion of roughness, a config¬ 
urational average over the quenched interface disorder is 
performed using an assembly of N p = 50 non-interacting 
particles. 


III. NUMERICAL RESULTS 

The origin of the exchange bias effect in FM/AFM 
nanostructures is attributed to the existence of uncom¬ 
pensated AF spins at the FM-AF interface^. For a 
spherical nanoparticle the interface is highly uncompen¬ 
sated due to its curved shape. However, if a cubical par¬ 
ticle is formed with faces parallel to the family {100} 
of crystallographic planes of a SC lattice, the core-shell 
interface is compensated (with exceptions at the corner 
sites), leading to a vanishing small exchange bias effect. 
Uncompensated spins in a cubical core-shell particle ex¬ 
ist provided the cube sides coincide with lower symmetry 
planes of the SC lattice. To achieve this we rotate the 


underlying SC lattice by an arbitrary angle </> about the 
[001] axis. The dependence of H c and H e b on the ro¬ 
tation angle </> is shown in Fig[2] Notice that for all </> 
angles, shown in Fig[21 the applied field remains paral¬ 
lel to the easy axis of the nanoparticle. The coercivity 



<t>(deg) 



FIG. 2: (Color on line). Dependence of low temperature (T = 
10~ 3 ) coercivity (upper part) and exchange bias field (lower 
part) of cubical and spherical particles on the rotation angle 
(</>) of the underlying SC lattice about the crystallographic 
axis [001]. 

and exchange bias field of the cubical particle depend 
substantially on the crystallographic orientation of the 
underlying lattice, while the corresponding fields of the 
spherical particle exhibit the expected rotational invari¬ 
ance (Fig@. 

Cubes with faces along {100} planes (</> = 0) have 
nearly twice the coercivity of spheres of the same size. 
The mechanism leading to the enhanced coercivity of the 
FM core- AF shell nanoparticles is the ’drag’ of shell- 
interface spins by the core-interface spins during magne¬ 
tization reversal . 10 ’ 11 . Consequently, the variation of H c 
follows the variation in the number of dragged AF spins 
at the shell-interface ( N s hif ), shown in Table Q] The re¬ 
duction of coercivity with rotation angle, as seen for C58, 
is a finite size effect. Namely, the dragging mechanism 
becomes less efficient as the fraction of dragged AF spins 
relative to the FM spins decreases. This situation occurs 
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when the 0-angle increases (N s hif/[N C + N ci f] = 0.651 
for 0 = 0° and 0.516 for 0 = 45°). 

On the other hand, the exchange bias effect of the cubi¬ 
cal particle increases (in absolute value) significantly for 
0 ss 45°, approaching the value of H eb for a sphere with 
the same geometrical size. This behavior is understood 
by the increase of the number of uncompensated spins 
of the rotated cube, which becomes comparable to the 
corresponding number for a sphere, as shown in Table 
HI Overall, the data in Fig [2] demonstrate that for cubi¬ 
cal and spherical particles with almost the same size (see 
Table Q]) , the same shell thickness and ideal FM-AF in¬ 
terfaces, cubical particles exhibit higher H c values and 
lower H eb values than their spherical counterparts. 

We consider next nanoparticles with roughness at 
the core-shell interface. In Fig(3] we show the low- 
temperature (T = 10 -3 ) hysteresis loops of cubes and 
spheres with rough interfaces. The underlying SC lat¬ 
tice is not rotated (0 = 0) in these calculations, thus the 
loops depict the effect of roughness on the compensated 
interface of a cube (C 58) and the uncompensated inter¬ 
face of the spheres (558, 5710). In all cases, roughness 
causes loop shearing, due to the development of a wide 
barrier distribution at the interface region arising from 
the mixing of the two phases. Notice that the shearing 
effect occurs for the average loop, which represents the 
mean behavior of an assembly of particles with different 
realizations of interface disorder. Individual loops (not 
shown here) have an almost rectangular shape, but they 
show a dispersion of H c and H eb values, which leads upon 
averaging to the observed shearing effect of the average 
loop.— Comparison of loops for 558 and 5710 indicates 
that the shearing effect is stronger for the sphere with 
the smaller core, because of the increasing interface-to- 
volume ratio with decreasing core size. 

In Fig [I] we depict the dependence of H c on interface 
roughness. For both cubes and spheres, H c shows a 
monotonous decrease with increasing roughness owing to 
the mixing of FM-AF at the interface that makes less 
effective the dragging of the shell interface moments by 
the core-interface moments during the magnetization re¬ 
versal of the FM. A similar reduction of H c with con¬ 
trolled interface roughness has been recently observed 
experimentally—. Cubes display higher coercivity than 
spheres of the same size, Hf 58 > H^ 5S , even in the pres¬ 
ence of roughness. Interesting, however, is the dramatic 
effect that even weak roughness can have on the coer¬ 
civity. As seen in FigU the presence of weak roughness 
(a « 0.05a) causes a dramatic drop of the coercivity of 
cube (758 below the coercivity of sphere 5710. The sen¬ 
sitivity of H c of the compensated cubical interface ((758) 
to disorder, implies that the experimentally measured co¬ 
ercivity values for cubes and spheres of similar sizes are 
screened by the presence of disorder. Detailed informa¬ 
tion on their interface structure^ 3 would be extremely 
valuable before any conclusion correlating the coercivity 
values to the nanoparticle shape could be drawn.— ^ 

Consider next the effect of roughness on the exchange 
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FIG. 3: (Color on line). Low-temperature (T = 10 -3 ) hys¬ 
teresis loops of cubical and spherical particles with rough in¬ 
terfaces. Error bars arise from configurational averaging over 
the interface disorder. 

bias field, as depicted in Fig(5] Interface roughness can 
either enhance or suppress the values of the exchange 
bias field. Particle (758 shows increasing values of \H eb \ 
with roughness, while (758r and the spherical particles 
show decreasing values of \H eb \. The distinct effect of 
roughness on cubes and spheres is related to the com¬ 
pensation of their interfaces. In particular, the increase 
of H eb with roughness in (758 is in accordance to the 
predictions of the random field model^, namely that 
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a/a 

FIG. 4: (Color on line). Dependence of coercivity on interface 
roughness at low temperature (T = 10~ 3 ). 
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FIG. 5: (Color on line). Dependence of exchange bias field 
on interface roughness at low temperature (T = 10 -3 ). 


roughness on a flat compensated AF interface, creates 
uncompensated spins required for the loop to shift. Con¬ 
trary to this trend, atomic scale roughness on uncompen¬ 
sated interfaces reduces the number of uncompensated 
spins leading to suppression of H e b , as observed on the 
flat interface of C58r and the curved interfaces of both 
spheres. As the degree of roughness increases, the shape- 
dependent Heb, is gradually suppressed and cubes and 
spheres of similar size have very similar values of H e {, 
for a > 0.3. Furthermore, an interesting feature seen in 
Fig© is the pronounced minimum of H e b at a ~ 0.05 
for C58r and S5 8. This stems form the wide dispersion 
of H e b values due to disorder that often assume posi¬ 
tive values and produce on average a minimum in |H e b|- 
The non-monotonous dependence of H e b on roughness 
has been previously observed^ - in exchanged-coupled bi¬ 
layers and interpreted as the competition between inter¬ 
face exchange and Zeeman energy of the AF spins at the 
shell-interface, when the interface exchange is antiferro¬ 


magnetic. We show here that a similar mechanism acts 
in core-shell nanoparticles irrespectively of their shape. 

The origin of the exchange bias effect in our model is 
the net (uncompensated) AF moment ( Maf ) that is ex¬ 
changed coupled to the FM^ Due to large anisotropy of 
the AF, which nearly freezes the net moment, the result¬ 
ing exchange bias field acting on the FM is expected to 
obey a linear scaling law, namely H e & = IMaf + H$. 
When an assembly of nanoparticle is considered, Maf is 
expected to assume different values due to different real¬ 
izations of the interface disorder. In this case, the scaling 
coefficient can be written as b = dH^IAM af , where the 
denominator expresses the variation of Maf values due 
to interface randomness. Therefore, one could interpret 
b as the average (over a nanoparticle assembly) sensitiv¬ 
ity of H e b to the microstructural details of the interface. 
In Fig. ©we depict the values of H e b for an assembly of 
non-interacting particles as a function of the net moment 
of the AF shell at the FC state. The linear fit to the data 
provides the average sensitivity parameter (6), which is 
found to increase as the the uncompensation of the AF in¬ 
terface increases. It is notable also, that larger dispersion 
of H e b values occurs for cubes than for spheres, because 
interface structural imperfections constitute a stronger 
perturbation to compensated (cubic) interface, namely 
an assembly of cubes is characterized by larger disper¬ 
sion of H e (,values than an assembly of spheres with the 
same size and degree of interface roughness. 



M af 


FIG. 6: (Color on line). Scaling behavior of exchange bias 
field with the net magnetic moment of the AF shell for spher¬ 
ical and cubical nanoparticles with rough interfaces. Data 
in each case are selected from an assembly of A/s = 150-200 
non-interacting nanoparticles with interface roughness in the 
range a =0.05-0.20. The slopes of the least squares fits, 
shown by straight lines, are b(S58) = 0.039(±2%), 6(6758) = 
0.008(±11%) and 6(C58r) = 0.017(±6%) 
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IV. CONCLUSIONS 

In conclusion, our simulations of the exchange bias ef¬ 
fect in FM core - AF shell nanoparticles with cubical 
and spherical shape show that among particles with ideal 
(non-disordered) core-shell interfaces cubes have higher 
coercive field (H£ ube > H^ phere ), owing to the larger 
number of ’dragged’ AF interface spins. This trend re¬ 
mains for uncompensated interfaces even in the presence 
of interface roughness. For compensated interfaces how¬ 
ever, weak disorder can reverse this situation rendering 
spheres harder than cubes. On the other hand, spheres 
exhibit higher exchange bias field (^jj^P here > H^ be ) due 
to the larger number of uncompensated spins on curved 
FM-AF interfaces. A non-monotonous dependence of 
H e h on roughness is observed, for small nanoparticles 
( R c < 5a) with antiferromagnetic interface exchange cou¬ 
pling and uncompensated interfaces, similar to the case 
of exchange coupled layered systems. Experiments on 
core-shell bi-magnetic nanoparticles with controlled in¬ 


terface roughness^! are required to investigate further 
this point. Strong roughness (a > 0.3) suppresses the 
shape-dependence of H e b and H c . The atomistic details 
of the interface are washed out. The H e b nearly van¬ 
ishes for both cubes and spheres and H c retains a weak 
shape-dependence. Shape-effects persist in the scaling 
factor relating the average exchange bias field of a (non¬ 
interacting) nanoparticle assembly to the net AF mo¬ 
ment, characterized by higher values for the spherical 
particles (b s P here > b cube ). Assemblies of cubical parti¬ 
cles are expected to have by a larger dispersion of H e b 
values than assemblies of spherical particles. 
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